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Abstract 

We compute the next-to-next-to-leading order (NNLO) soft and virtual QCD 
corrections for the partonic cross section of colourless-final state processes in hadronic 
collisions. The results are valid to all orders in the dimensional regularization 
parameter e. The dependence of the results on a particular process is given through 
finite contributions to the one and two-loop amplitudes. To evaluate the accuracy of 
the soft-virtual approximation we compare it with the full NNLO result for Drell-Yan 
and Higgs boson production via gluon fusion. We also provide a universal expression 
for the hard coefficient needed to perform threshold resummation up to next-to-next- 
to-leading logarithmic (NNLL) accuracy. 
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1 Introduction 



The development of accurate QCD calculations is a fundamental tool to properly test the Standard 
Model. Given the size of the perturbative corrections, leading order (LO) evaluations are insuffi- 
cient, and higher perturbative orders must be taken into account. However, these calculations are 
highly non-trivial, and at present only a few processes have been computed analytically with full 
next-to-next-to-leading order (NNLO) precision. At hadronic colliders, only Drell-Yan [TJ [2] and 
Higgs boson production [21 El H] (within the effective vertex approach) have reached that stage of 
accuracy. 

Higgs production is a particular example of an observable with a slow convergence for the 
perturbative expansion in the strong coupling constant as- Next-to-leading (NLO) corrections 
[HI El [7] are as large as the Born result and the NNLO contribution still increases the cross section 
by about 25% at LHC energies. Analyses from scale variations [HI El Ell EH EG] and soft-gluon 
expansion 

P31 HH Q31 QUI 02] indicate that the next orders (N 3 LO and beyond) can contribute 
still at the level of 10%. 

An important step towards a complete NNLO calculation for both Drell-Yan and Higgs pro- 
duction has been the evaluation of the soft and virtual contributions [T8 l fT9 l 120] . which provide the 
dominant terms for those processes. As a matter of fact, this is a general feature when a system of 
large invariant mass Q is produced in hadronic collisions. Since parton distributions f a /h( x ) grow 
very rapidly for small fractions of the hadron momentum x, the partonic center-of-mass energy 
tends to be close to the invariant mass Q, and the remaining energy only allows for the emission 
of soft particles. For this reason the soft-virtual contributions are expected to be a very good 
approximation to the total cross section for a large number of processes. 

In this paper we exploit the factorization properties of the QCD matrix elements to compute 
the soft and virtual contributions to the partonic cross sections at NNLO for a wide number 
of processes in hadronic collisions where a system of colourless particles is produced (as gauge 
bosons, Higgs, leptons, etc.). The computational approach presented here can be extended to 
higher orders in perturbation theory simplifying considerably the evaluation of the corrections. 

We present a universal expression for the corresponding cross section at NNLO, valid to all 
orders in the dimensional regularization parameter e. With this result, it is possible to evaluate 
the soft- virtual approximation for any process of the kind studied in this paper in an automatized 
way once the relevant one- and two-loop amplitudes become available and, therefore, provide a 
first estimate of the size of higher order corrections for a number of interesting observables at the 
LHC. Furthermore, since the results are valid to all orders in e, they become necessary ingredients 
for ultraviolet and infrared factorization at N 3 LO (and beyond) within the same approximation. 

Another interesting use of the soft approximation is the relation to the soft-gluon threshold 
resummation approach. We profit from this calculation and obtain, for the first time, a universal 
expression for the hard coefficient needed to perform threshold resummation up to next-to-next- 
to-leading logarithmic (NNLL) accuracy. 

The paper is organized as follows. In section [2] we present the notation for the QCD cross 
sections and show how phase space factorization occurs in the soft limit. In sections [3] and H] 
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we perform the calculation of the soft-virtual corrections at NLO and NNLO respectively. In 
section [5] we present the soft-virtual approximation in Mellin space. In section [6] we analyse the 
phenomenological results for Drell-Yan and Higgs boson production via gluon fusion to compare 
the soft-virtual approximation with the full result as a way to validate its dominance. In section 
[7] we profit from the previous result and present a universal expression for the hard coefficient 
required to perform threshold resummation up to NNLL accuracy. Finally, in section [8] we present 
our conclusions. 



2 QCD cross sections 



We consider the following general process in hadronic collisions: 

h + hz^F + X (1) 

where F denotes any colourless final state (i.e. without quarks or gluons), and X stands for any 
inclusive final hadronic state. The center-of-mass energy is y/sjj, and Q is the invariant mass of 
the system F which can involve a combination of gauge bosons, Higgs, isolated photons, leptons, 
etc. The inclusive cross section can be written as 



do ( 
Q 2 Ira^^) = S / d%l d%2 f*/hi( x i> fb/h 2 (x 2 , fi 2 F ) / dzb[z- 
"v „ , Jo Jo V 



x±x 2 



a,b 

xa z G ab (z; a s (/4), Q 2 /v R ] Q 2 7 74) , ( 2 ) 

where r = Q 2 /sh, and hr are the factorization and renormalization scales respectively, and &o 
is the Born level partonic cross section. The parton densities of the colliding hadrons are denoted 
by fa/h( x , A*f) an d the subscripts a, b label the type of massless partons (a, b = g, qf, g/, with Nf 
different flavours of light quarks) . 

According to Eq. (j5J), the cross section a a b for the partonic subprocess ab — > F + X at the 
center-of-mass energy s = X1X2SH is 

Q 2 ^(s,Q 2 ) = - s a Q 2 G ab {z) = a z G ab (z) , (3) 

where the term 1/s arises from the flux factor and leads to an overall z factor, being z = Q 2 /s 
the partonic equivalent of r. The hard coefficient function G ab has a perturbative expansion in 
terms of powers of the QCD renormalized coupling as{fi 2 R )'- 

G ab (z;a s ^ 2 R ),Q 2 /^Q 2 /ti) = ^ ^ ■ (4) 

n=0 ^ ' 

In the following, the dependence of «s on the renormalization scale is understood. We always 
use the MS scheme for the renormalization of the strong coupling. 

At leading-order the partonic subprocess is ab — > F, and since the final state is colourless we 
can only have ab = gg or ab = qq (and ab = qq). The LO contribution is then 

GabW = S a 9 h 9 6(1 - z) oi G®(z) = (6 aq 5 bq + 5 a , 5 bq ) 6(1 - z) . (5) 
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At higher orders, other parton subprocesses can contribute to the total cross section. 



In the soft-virtual approximation, however, we are only interested in the same parton sub- 
process present at LO, and we compute only those contributions in G a a that give rise to the 
distributions 5(1 — z) and T>i(z) in the coefficient function, where we have defined 



The + indicates the usual plus-prescription, 



ln*(l - z) 



(6) 



1 f 1 lnYl - z) 
V i (z)f(z)dz = — } ) -[f(z)-f(0)]dz. 

JO 
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(7) 



These two types of contributions are the most singular terms when z — > 1, and then dominate the 
cross section in the soft limit. 



The aim of this paper is to calculate the NNLO soft-virtual corrections for any process of the 
type of Eq. (Op). The parton subprocesses that contribute up to second order in as are 





aa — > F 


(tree-level) 




aa — » F 


(one-loop) 




ad — > F + g 


(tree-level) 


q| 


aa — > F 


(two-loop) 




aa — >■ F + g 


(one-loop) 




aa — >■ F + q + q 


(tree-level) 




aa — >■ F + g + g 


(tree-level) , 


with aa = gg or aa = qq, depending 


on the process. 




2.1 Phase-space factorization 




To compute the real corrections of the processes we 


are interested in 



(8) 



phase-space integration of the corresponding matrix elements in the limit in which the emitted 
QCD-partons become soft. In this limit, the phase-space can be written in a factorized form. 

Let us consider that the final state F has / non-QCD particles with momenta {k}, and m soft 
QCD massless partons with momenta {q}. The n-dimensional phase-space is then 



/^/ 



m Jn 



h + . . . h E 



( y 2n) n 5 n (p 1 +p 2 -qi- . . . q m -ki- . . . h 

(9) 



We introduce the momentum K 
by the identity 

1 = J d n K 5 n {K -kx-...h) J dQ 2 5 + {K 2 - Q 



with K = Q . Multiplying the above equation 



(10) 
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we arrive at 



dPS 



dQ 2 
2tt 



.i=i 



{2i:) n 8 n {K - h - . . . k) 



x 



m A n 



d n K 



L(2tt)' 



-5+(K 2 -Q' 



{27T) n 8 n {p 1 +p 2 -K-q 1 



q m ) ■ 



In the soft limit we have 



5 n (K-h- ...h) 5 n (p 1 +p 2 ~K-q 1 - . . . q m ) ~ 8 n {p 1 +p 2 -k 1 - ...hi) 8 n ( Pl +p 2 -K- qi - . . . q m ) , 

(12) 

g 

where the symbol ~ indicates that the equality is valid when the emitted QCD partons are soft. 
Within this approximation we obtain in the first line of Eq. (llljl the corresponding leading-order 
phase-space c?PS^ \ which contains the dependence on the internal variables of the system F. The 
second line in Eq. fTTT]) is the phase-space of a process with one particle of invariant mass Q in 
the final state plus m soft partons, dPS 2 ^ 1+msoft Then Eq. (TH]) can be rewritten in the following 
way: 

dQ2 f dPS® I dPS 2 ^ 1+mso{t , (13) 



rfPS 



2ti 



arriving to a factorized expression for the phase-space in the soft limit. 



3 NLO 



At NLO we have to consider the one-loop corrections to the partonic subprocess aa — > F, and 
also the real gluon emission subprocess, ad — > F + g. We begin by computing the latter. 

Let Ai^ be the LO matrix element, and Aif^ the correction corresponding to the real gluon 
emission subprocess at tree-level. In the limit where the momentum q of the gluon becomes soft, 
M-f^ can be written in the following factorized way |22j : 

\Mf\q, Pl ,p 2 )\ 2 ~ (ao^)8rrS g (q,p u p 2 )C a \M^( Pl ,p 2 )\\ (14) 

where p\ and p 2 are the momenta of the incoming QCD-partons, and the dependence of the matrix 
elements on other non-QCD particles momenta is understood. The symbol «o stands for the bare 
coupling constant, and Hq is the dimensional-regularization scale. Renormalization is achieved by 
the replacement 

«o vl £ S e = «s /4 e - «s y + ^("s)^ ; ( 15 ) 

where (3 is the first coefficient of the QCD beta function and S e is the typical phase-space volume 
factor in n = 4 — 2e dimensions: 

A) = UCa ~ 2Nf , S e = M'e-r*, (16) 

1Z7T 

being ^ E = 0.5772 . . . the Euler number. 
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The scalar eikonal function S g (q,Pi,P2) contains all the dependence of M.^ on the soft gluon 
momentum, and takes the form [22] 

Sg{q,Pi,P2) = 7 w r, (17) 

while the constant C a depends on the nature of the radiating parton, being 

C a 



C A = N c , if a = g 

C F = (N 2 c -l)/(2N c ), if a = q,q 



where N r = 3 is the number of colours. 



Combining Eq. ([14]) with the phase-space factorization of Eq. ffTBj) we arrive at the following 
expression for the NLO tree- level real gluon emission cross section a g ^: 

j-(o) ~ r 

aa g s (T / ^ / c , \ tt3q2 — ^l+i soft 



dQ 2 2vr 



(a ^)MC a J S g (q, Pl ,p 2 )dPS 2 ^ 1+lsoit . (19) 



The phase-space integration in Eq. (fT9j) can be performed in a closed form. After some simple 
algebra we arrive at 

Using polar coordinates in n dimensions we can write 

d n ~ l q = d\q\ \q\ n - 2 dn n _ 2 , (21) 

with 

dfi n _ 2 = sin"" 3 0! sin"" 4 2 . . . sin fl„_ 3 d0 2 . . . d0„_ 2 , (22) 

where n _ 2 G [0, 2it) and 0j G [0, 7r] for the others. Parametrizing the momenta in the center-of- 
mass of the incoming partons, and setting n = 4 — 2e, it can be shown that 

S g (q, Pl ,p 2 )dPS^°* = — S^g / e d9 , (23) 

s i+£ (l — z) l+Ze 27rr(l — e) J (1 — cos#)(l + cos#) 

where 9 = 0i is the angle between the soft gluon and the n-th axis. Using the variable y = 
(1 + cos6 l )/2, the remaining integral can be carried out as a particular case of 

o l(a + pj 

Finally, the tree-level real gluon emission cross section in the soft limit has the following 
expression: 

ddf s, _ao /Op\ (47rjzg/ 8 ) e 4T(l-e) 

dQ 2 " 8 \2 J (1 - a e T(l - 2e) ' K ' 

This formula is valid to all orders in e for any reaction of the kind of Eq. ([I]), and its only 
dependence on a particular process is in the Born- level cross section a . The C a factor only 
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depends on the nature of the incoming partons. The expansion of the factor (1 - 
the appearance of 5(1 — z) and T) i (z) terms, according to the following relation: 

1 1 . v^f-aeV' 



-1-2<E 



1 - z) 1+ae 



ae 



5(1 



i=0 



V % (z) . 



leads to 



(26) 



We now have to evaluate the one-loop correction term. Even though this cannot be done 
in a process independent way, the infrared-singular behaviour of QCD amplitudes at one-loop 
(and two- loop) order is well known [23l EH [25J [26]. For the processes we are interested in, the 
renormalized one- loop order amplitude Ai^ can be written in terms of the Born- level amplitude 
M {0) in the following wajQ [26]: 



where 



2vr 



S7 1 



„ 1 1 
C n - + 7a - 
e 



e l 



(27) 



(2f 



s j r(i 

and the contribution Aig^ is finite when e — > 0. The coefficient 7 a depends on the initial-state 
partons, being 



7« 



~e A 
\C F , 



\N f , \ia = g 



if a 



q,Q 



(29) 



For the NLO calculation we only need the 0(a^) term of the squared matrix element, that is 
Ai^ (Ai^)* + (Ai^)*A4^. Performing the formal phase-space integration of this term we arrive 
at the following expression for the one-loop virtual contribution to the cross section: 



da 



(l) 



S7 i m f 1 i 

r(l-2e)r(2e) V a e 2 e 



0o 



dQ 2 

where (e) is a one-loop finite contribution to the cross section defined by 



O"0 



(30) 



(31) 



To obtain the NLO soft- virtual contribution to the coefficient function G a a{z) we still have to 
add to Eqs. (1251) and (|30|) the counterterms coming from mass factorization. Keeping terms up to 
second order in powers of as, relevant for the NNLO calculation in the next section, we have 



da(s) 



s I 2ir 



G^(z)--P£\z) 



+ 



QtsV 

2tJ 



2K ) ®Pi° ) )W+2 7 r/3 Pi 0) W 



(32) 



p£>(z) + 2(p<»®g. 



+ G%>(z)} + (D(a s 



TFor the sake of simplicity we omit the explicit dependence of the matrix elements on the partons momenta. 

"'"We explicitly keep higher order terms in e originated in the one-loop amplitude as they contribute to the final 
result to the same order in the dimensional regularization parameter. Those in are implicitly included in the 

definition of <t q = ^ / {M^ | 2 dPS (0) . 
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where for simplicity we have set fip = Q. The symbol "£g>" stands for the usual convolution. The 
Altarelli-Parisi splitting functions in the soft limit (z — > 1) take the form (27] 



P®(z) ~ la S(l-z) + 2C a V (z 
P$t(z) * ^5(l-z)+C a 



... . 67 tt 2 \ Ar 10 



where we have defined 

C| (| + 3C 3 ) - ±CViV> - fCUty, if a = g 

Cl{l-^ + ^) + C F C A {^+ l ^-Kz)-C F N f {^ 2 + ^), ifa = g,g 



la 



(33) 
(34) 



(35) 



Combining the results of Eqs. ( 125]) and ( 1301 with the counterterms coming from Eq. ( |32i) we 
arrive at a closed expression valid to all orders in e for the NLO soft- virtual coefficient function 
G a ^. For simplicity, we only write the first three terms of its expansion in powers of S 



G2M ~ 5(1 -z) Jc.?£ + Sa + 

3 (T 



+ 



-(1) 

- 4G a C3 + 7a— H ; — 

6 a 

.(1) 



r 177r \, 2 C 3 , ff M2) 
C„— tt— + 7a^T H : 



90 



00 



(36) 



28Cs 



16. 



+ 8C fl D 1 ( Z ) + C a ^7^0(2) - 8£> 2 (z)Je + C a ^-^D (^) - 2^1(2) + —V z (z)y 2 + 0(e d ) 

s 

where we have set lir = lif = Q — y/s. The dependence of this expression on a particular process 
is contained only in the one-loop coefficients (Tg^(i) defined by 

00 

^) = £4V- ( 37 ) 



8=0 



4 NNLO 



At second order in the perturbative expansion we have to consider the four 0(a^) parton subpro- 
cesses of Eq. ([8]). The one-loop correction to the subprocess aa F + g can be obtained in a very 
similar way to the tree-level one. Let M.^ be the one-loop gluon emission matrix element. The 
analogous to Eq. (fT4l) is given by the soft limit of one-loop amplitudes as [281 ISH EUJ (3TJ [32] 

(M {1) )* M (0) + M (1) (M (0) )*~ 

- c * ¥ 2 ^|i'5ff" ^"» I+ ^ ,0 T} ■ < 38 > 

''Contributions up to C(e 2 ) are needed to build the renormalization and factorization counterterms for a calcu- 
lation to N 3 LO accuracy. 
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The phase-space integrals we have to perform are 

S g (q,p uP2 )dPS 2 ^ 1+lsoit and / [S g (q, Pl ,p 2 )] 1+e dPS 2 ^ 1+lsoft 



and their calculation can be achieved with the tools discussed in section [3j The result is 



da. 



(i) 



dQ* 



S 



/a \ (47r/xg/s) e _ 4T(l-e 



s \2tt7 (l-z)U- 2e a e r(l-2e) 

e)r 3 (e)r( 2,) 



(39) 



(40) 



_ ao /aoy (fTT^/sT r 2T 3 (1 

s V2J a-z)^ a A rri- 



4e)T(l - 2e)r(2e) 



where a^ can be written as in Eq. ( )30l) . using the delta function to perform the Q 2 integral. 

We continue by computing the double real emission subprocesses, that is act — >■ F + g + g and 
aa —7- F + g + p. For the NNLO squared amplitudes the following infrared factorization formulae 
hold 



Mqq(q 1 ,q2,Pi,P2)\ 2 ^ ("0^) * ^ C * [ 2*1 (gl > 92) + Z22 (gl , ^2) - 2Xi2 (gi , g 2 ) ] I A4 (Pi , P2) 1 2 , (41) 



,2e\ 2 o„2 



(42) 



A4 9g (gi,g 2 ,pi,p 2 )| ^ ("oA^) 16tt 2 |.Mo(Pi,P2)| 

AClS g (q{) S g (q 2 ) - C A C a (Sii{qi,q 2 ) + S 22 (gi,g2) - 2 <Si 2 (gi,g 2 ) 

where all the dependence on the momenta gi and q 2 of the soft particles is embodied in the 
functions and 5^-, which take the form [33] 



2y(gi,g 2 ) 
^j(gi,g 2 ) 



+ 



2(p» ■ Pj) (gi • 92) + [p» • (gi - 92)] [Pj • (gi - 92)] 
2(gi • g 2 ) 2 • (gi + g 2 )] \pj ■ (g x + g 2 )] 

'1 - e) pj • qx pj ■ q 2 +Pi-q 2 pj ■ gi 



(?i • g2) 2 Pi • (gi + 52) Pj • (<?i + 92) 



(p< • Pi) 



2pi • gi Pj • g 2 Pi ■ g 2 Pj • gi 



2 - 



Pi • qiPj ■ g 2 + Pi • g 2 Pj • gi 

Pi • (gl + g2) Pj • (gl + Q2) 



Pi ■ Pj 



+ 



2gl ■ g 2 LPi • gl Pj ■ g2 Pj • gl Pi • g2 



Pi ■ (gl + g2) Pj • (gl + g2) 
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(gj • glPj • g2 +Pi • g2Pj • gl) 
Pi • gl Pj ■ g2 Pi • g2 Pj ■ gl 



(43) 



(44) 



(45) 



To perform the phase-space integration of both contributions we use the parametrization of 
^pg2->i+2soft introduced i n [18] for the calculation of the second order corrections to the Drell-Yan 
process. Introducing the variables g = gi + g 2 and si 2 = g 2 we can write 



dPS 



2->l+2 soft 



(2tt) 2 "- 3 

x J d n K f d n q [ ds 12 5 + {K 2 - Q 2 )5 + {q 2 - s 12 )5 n { Pl + p 2 - K - g) 

x 



J d n Ql J d n q 2 o+{ql)5 + {ql)5 n {q- qi -q 2 ) . 



(46) 
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The last line of the above equation is most easily computed in the center-of-mass of q% and q 2 . In 
this frame, orientated so that p\ is in the direction of the n-th axis and p 2 lies in the plane defined 
by the n-th and (n — l)-th axes, the momenta can be parametrized as follows [18J: 

Qi — |\/si2 (1j • • • ^ cos sin 6 1 , cos 9) , 
q2 = |a/^i2 (1, • • • , — cos sin 6 1 , — cos#) , 
(s-t) , 



' s — Q 2 _ Sl2 

if = ( — — , 0, . . . , 0, 1 If I sin^>, \K\ cosip 



A/A(g,Q 2 ,gl 2 ) 

2J^~ 2 



|if| 

co = (.-<y)( 6 -g')-. a (f + g') 

l>-f)VAl>,QVs 12 ) 
where we have defined 9 = 9\ and (f) = 9 2 (see Eq. f[2"2"j) ). and 

t = 2 Pl -K, u = 2p 2 -K, s = (p 1 +p 2 ) 2 , s 12 = s-t-u + Q 2 , (48) 

where A (a, b, c) is the Kallen function, A (a, b, c) = a 2 + b 2 + c 2 — 2ab — 2bc — 2ca. The dots in 
qi and q 2 represent n — 3 unspecified components of momentum, which are trivially integrated 
with the present parametrization. Using momentum conservation, we obtain for p 2 the following 
expression: 

S\ 2 ~\~ t Q 2 , V 

p 2 = — — — (1,0, . . . ,0,sinx,cosx) , 



cosx = — - ( |if|cos^ + ^ ) . (49) 

s l2 + t-Q 2 V 1 ^ 2y/S^J 



Considering the above parametrization and introducing the variables x, y and z: 

z = Q 2 /s, 
u = s(l-y(l-z)), 

( , s 2/(1 — 3/)a;(l - z) 2 \ , . 

the double real emission phase-space can be written as follows |18j : 

i „n— 3 pit pit 

dPS 2 ^ 1+2soft = -V — (1 - z) 2n ~ 5 d9 d(j) sm n ~ 3 9 sin n - 4 

(4tt)« r(n - 3) 1 J Jo Jo 

X / / [y(l - y)]- 3 [z(l - x)]™/ 2 - 2 [1 - y(l - z)] 1 ^ . (51) 



o Jo 



Up to this point we have kept the exact expression for the phase-space. In the soft limit the 
last factor of Eq. (IBTjl can be approximated by 1, obtaining 

d pg2^1+2soft £ ^ 3 Q _ z )2n-5 H H - £ sin «-4 - 



(47r)»r(n-3) Jo 
x rfx [y(l - y)]- 3 [ar(l - x)]"/ 2 - 2 . (52) 



JO 



9 



With this approximation the integrals we have to perform become considerably simpler. 

We now have to express the integrands in terms of the variables x, y, z, 9 and 0. The quark- 
antiquark emission subprocess is relatively simple. The Xn (51,^2) term of Eq. (HI]) , for example, 
takes the following form 

2 cos 2 6 _ 2(1 - (1 - z)y) 2 cos 2 6 ^ 2 cos 2 9 

2-n Ui, 92) - ^-y- - -^jiT^Y[iqjv ~ ~ s 2 (l- z) i y 2 (l-y) 2 x 2 ' ^ ' 

where we have considered the z — > 1 limit in the numerator in the last step. Then the phase-space 
integral can be done in a direct way, obtaining 

X (a fl UPS 2 * 2 »^ r(-l-e)r(-2e)r(2-e) 
X n ( gi ,g 2 )rfPS -- a i + 2 e(1 _ ^ 32vr3r(-4e)r(4-2 e ) " (54) 

The phase-space integral of the X 22 (gi, 92) term yields the same result, since the integrand can be 
obtained from X n by the exchange pi ++ p%. The same happens for two integrands that only differ 
in the exchange of q\ ++ q<i. 

The remaining term is Xi2(gi, 92)- It can be split into two contributions, X12 = X4 +Xb, where 

X A (q x q 2 ) = - Pi) (gi • gg) £ 4 

X ' 2 (91 ■ 92) 2 [pi • (91 + ^2)] [P2 • (91 + 92)] s 2 (1 - z) 4 y 2 (l - y) 2 x ' 

2^(9 9 ) ~ ~ f pi ' ( gl ~ g2 )J ^ 2 ' ( gl ~ g2 )] ~ 2x cos t sin 26> ~ ( 2 ~ 4x ) cos2 ® ^ 56 ) 
2(9i • 92) 2 bi • (91 + 92)] [j?2 • (91 + 92)] s 2 (1 - z) 4 y 2 (l - y) 2 x 2 

Their integration is straightforward and the result is 

XAa a )dPS™+ M ~ {4e)2e r(-e)r(-2e)r(l-e) 

X A ( gi ,g 2 )rfPS - - a i+2e (1 _ ^1+4, l 67 r3r(l-4 e )r(2-2e) ' (57) 

n „ fl W p S 2^i + 2sof t s (4tt) 2 - r(-l-e)r(-2e)r(2-e) 

X^^PS _ - s i +2e(1 _ 3)1+4, 87 r3r(l-4 e )r(4-2e) ' (58) 

Combining Eqs. ( |54"1) . ( 1571) and (1351) with Eq. (T4Ti) we arrive at a rather compact expression for 
the NNLO double real soft quark-antiquark emission cross section: 

dta b, (ao\ 2 (47r^/gf 2 r 2 (2- e )r(-26) 



dQ 2 " s V27r7 (l-z) 14 - 4e C ' aiV/ e 2 r(4-2e)r(-4e) ' (59) 
The Nf factor arises from the sum over the different qq pair flavours. 



The integration of the double gluon emission matrix element of Eq. (1421) is more cumbersome. 
The following formula is used several times in the calculation of the angular integrals [18J: 



, sin"" 3 9 sin"" 4 6 , . 

' (1 - cos6») ? (l - cos£cos6» - sm^sm^ v ' 







2 i-,- % r(n/2-l-j)r(n/2-l-») T(n - 3) / n f 
T(n-2-i-j) r 2 (n/2-l) V' J '2 ' S 2, 
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where F(a, b; c; z) is the hypergeometric function. 

As an example we show in detail the integration of the simplest term, J S g (qx) iS 9 (g 2 )^PS 2 ^ 1+2soft . 
After partial fractioning, it can be written in the following way 

S g { qi )S g {q 2 ) = - S \ + + . (61) 

(s -t) (sia +t- Q 2 ) \Pi ■ qi Pi-Q2j\P2-qi Vi-qi) 

Out of the four terms we obtain multiplying the factors in parenthesis we have to integrate only 
two because the others can be obtained from them by the exchange p± -H- p 2 or q± <H- q 2 . 

Let Q2) be the term corresponding to the product (pi • qi){j>2 • 9i) in the denominator. 

Using the parametrization of Eq. (1471) we arrive at 

Sa(Qi, ?a) = (( s _ l ) ( Sl2 + f _ g 2 ) ) ((1 _cos0)(l -cosxcos^-sinxcos^sin^)) ' ^ 

The first factor can be written in terms of the variables x, y and z. Then the second factor can 
be integrated using Eq. fl60l) with i — j — 1 and £ = x (the factors sin™ -3 sin n_4 cf> are present 
in c?PS 2 ^ 1+2soft ). The argument of the hypergeometric function, cos 2 |, can be written using the 
definition of cos \ in Eq. (|49|) in the following way: 

2 x 1 + cos x s , . 

cos — = ~ 1 + x , 63 

2 2' v ; 

where the last step is valid in the soft limit, since cosx — > 1 — 2x when z — > 1. 

Now we have to perform the integration in the variables x and y. The first one takes the 
following form in the soft limit 



I F(l,l;^-l]l-x)[x{l-x)] n / 2 - 2 dx 
Jo 



which is a particular case of 

r 2 ( 7 )r(2 7 - 2) 



/ f(i,i, 7 ,*)[(i- t)ty- x dt 

Jo 



(64) 



(65) 



r 2 (2 7 -i) ' 

The remaining integral can be done straightforwardly, and the final result is 

W*.*)***™-*" aJTrpU) • (66) 

Now let «Sfi(gi, 92) be the term of Eq. ( 1611 corresponding to the product (pi • qi){j>2 • 92) in the 
denominator. In this case we have 

& (gi, ga) = ( ( s _ £ ) ( Sl2 + l_ g 2 ) ) ((1 -cos#)(l + cosxcosfl + sinxcos0sin#)) ' ^ 



To perform the angular integral we use again Eq. ( 1601) . taking in this case i = j = 1 and £ = X + 77 - 
Then the result only differs from the previous one in the argument of the hypergeometric function, 
which is 



2 £ 2 I X + k\ 1 - cos* * 



cos - = cos ^ — - — ) = ~ x , (68) 
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and the x integral takes the form 



F(l, 1; f - 1; x)[x(l - x)] n/2 - 2 rfa; , 



(69) 



which is the same of Eq. (1641) through the change of variables x — > 1 — x. Then the phase-space 
integral of J <Sb(<Zi, 92) <iPS 2 ^ 1+2soft coincides with Eq. ( )66|) . and the four terms of Eq. ( )6T|) give 
equal contributions to the cross section, yielding 



S g ( qi )S g (q 2 )dP$ 



2-S-1+2 soft f, 



(47T 



1 2c 



r 2 f-el 



s i+2e(i _ z y+*e 87r 3 e r(l-4e) ' 



(70) 



The integration of the remaining terms of Eq. ( I42p is more complicated due to the more involved 
expression of the two-gluon eikonal function (Sy^i, 52)- However, the techniques needed for the 
calculation are essentially the ones we have described so far. The result we obtain is the following: 



J <Sn(<7i, q 2 ) + <S 2 2(?i, 92) - 2<Si 2 (gi, (h) d¥$ 



2^1+2 soft 



(4tt 



,2c 



r(-e)r(-2e)r(2 - e) 



s i+2e(i _ z y+*e I67r 3 e 2 r(-4e)r(4 - 2e) 



X 



- 3 + (19 - lle)e + 2(3 + 4e(e - 2)) 3 F 2 (1, 1, -e; 1 - 2e, 1 - e; 1) 



where the generalized hypergeometric function 3F2 arises from the following integral: 



(1 - xy'x-^Fil, 1; 1 - e, x)cix 



2T(-2e) 



} F 2 (l,l,-e;l-2e,l-e;l) . 



(71) 



(72) 



Combining the results of Eqs. (1701) and (ITT]) , adding the factors coming from Eq. (142]) and a ~ 
factor for identical particles in the final state, we finally arrive at the double soft gluon emission 
cross section: 



da, 



mi 



dQ* 



a /q \2 (47r/i 2 /g) 2e r(-e) / ^ 2tV ^ , ^ r 
TfeJ (l_ 2 )i + 4c e2r( _ 4e) |2C a r(-e) + C a C, 



r(-2e)r(2 - e) 



T(4 - 2e) 

3 + (19 - lle)e + 2(3 + 4e(e - 2)) 3 F 2 (1, 1, -e; 1 - 2e, 1 - e; 1) 



Again, this result is valid at all orders in e 0. 



(73) 



We now have to evaluate the 0(a§) virtual corrections. One contribution comes from the 
square of the matrix element Ai^ of Eq. ( 127|) . The other arises from the two-loop (renormalized) 
matrix element Ai^ 2 \ whose infrared-singular behaviour can be written as [26, 36J: 



Z7T \27T/ 



lV(e)MW(e)+M$(e) 



(74) 



^The contribution from the generalized hypergeometric function 3 F 2 is the only term that can not be written 
in terms of simpler T functions [35 . 
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where M.^ is finite when e — > 0, and the function I^\e) is explicitly given by [26] : 



72e 4 r(l-e) [ 



Ue(C a + ej a )(UC A -2N f )-36 



sr 1 



r l-e 



4vr/i 2 ^ e 



(C a + e 7a ) 2 (75) 



+ ef-1) 



/4 



36e 2 # a + 2(3 + 5e)(C a + 2e Ta )iV / + C A {C a + 2e 7a )(-33 - 67e + 3e7r 2 



Here the coefficient H a depends on the type of the incoming partons, being 



' C\ (|Ca + h + ^ ) - CUiV/ (f + g) + 1CWV) + | AT) , if a = <? 



C^(-6C 3 -l + f )+C A C F (fC3 + 



13 a i 245 23tt 2 



216 48 



+ C F N f 



24 108 / ' 



(76) 



if a = g, q 



The corresponding contribution to NNLO arises from the product with the Born-level matrix 
element, i.e., M {2) (M {0) )* + (M {2) )*M {0) . 

Combining all the second order virtual corrections we arrive at the following expression for the 

(2) 

two-loop virtual contribution to the cross section a\, . 



S- 1 /4ttm 2N ' 



72 e 4 



4vr/4 \ e WS-\C a + e la ) 2 T\e) 
T 2 (2e) r 2 (l -2e) 



+ e 



r(2e)r(l-2e) 



r(4e) T(l - e) T(l - 4e) 
x \C A (37r 2 e - 67e - 33) (C a + 2e 7a ) + 2N f (5e + 3)(C a + 2e 7a ) + 36H a t 

12(C a + e la )T(e] 



r(2e)T(l-2e) 



& {1) (e) 
llC A + 6e^^ + 2N f 



(TO 



(77) 



(2) 

where a^ n (e) is a second order finite contribution to the cross section defined by 

»2w . 



(TO 



A<£ (-M (0) )* + (MS) -M (0) ) + KnT dPS<°> / / |A*<°f dPS<°> . (78 



We have evaluated all the NNLO corrections for a process with a colourless final state. The 
results are valid at all orders in e. To obtain a finite result we must add Eqs. ( 1401) . ( 1591) . ( 1731) . 
f!77|) and the 0(a|) of Eq. (1251) together with the counterterms coming from mass factorization 
(Eq. fl32l) ). In this way we arrive at a closed expression for the second order coefficient function 
Gal (z). Expanding the result in powers of e and keeping terms up to (9(e) (relevant for an eventual 
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calculation at the next order) we obtain the following result: 

2tt 1 ... ... /607 , 737vr 2 407C 3 7tt 4 

~36 "48" 



*(1-*)<C£ — +C a C A (— + [32 



+ 



15 

12 55tt 2 37C; 



+ 



+ e 



C„ 2 



80tt 2 C: 



1 216 1 

3 



llTT 3 „ 27T 2 <5- 



(1) 



fin(O) ^(O) 



(2) 



O (JQ 



+ 



(TO 



+ CqCa 



s 7 ° f \ 54 243 81 960 J 
4087C 3 + 101vr 2 C 3 43C 5 + 3644 707tt 2 847tt 4 



108 



243 162 1920 



7vr 4 



a 



(i) 

fin(O) 



(1) 

fln(l) 



+ C a l_ 7a -4C 3 ^+ 

9 (X 3 cr 



+ C A la 



11C 3 1675vr 2 25tt 4 



; [a 1 9 108 



U 1 ) 



18 

(2) 



216 



72 



25tt 2 77r 2 7a OW (0) a fin(1) 



12 



6 ct 



(TO 



16tt^ 



/ 404 


+ : 


V 27 


/56 


An' 


V27 ~ 


9 


10tt 4 . , 
9 V Q {z) 



+ C 2 a [64( 3 V (z) - —T>i(z) + 32V 3 (z 

+ C a Ca 

+ Ca Nf 

< <cl 



268 4vr 2 
~9 3~ 



40 



V Q {z)-—V 1 {z) + -V 2 {z) 



a 



(i) 



8C a ^V 1 (z) 



(TO 



116tt 2 _ , , 160 _ , , 

V 2 (z) VAz 

3 y ' 3 V ' 



3vr 2 536 



1616 121tt' 



27 



9 



3 *) 



_ , . 88„ , . 

Z> 2 (*) + y© 3 (z) 



32C 3 328 35tt^ 
~3~ + ~81 27" 



V (z) + 



22vr 2 224 



+ Ca 



• <T (1) 

TT — L>o{Z) 

O"0 



287T 2 

3 cr 



a fin(l) A 



9 27 



_ , , 80_ , . 16_ , , 
Vi(z) + -V 2 (z) - -V 3 (z) 



V l {z) -8^^V 2 (z) 

O"0 



+ C(e 2 ) 



(79) 



where for simplicity we have set fj,p = = Q. The dependence of the above expression on a 
particular process is embodied in the finite virtual contributions o , ^r i \ and fr^u)-* the latter being 
defined by 



>=£■ 



fin(i) 



(80) 



i=0 



The coefficients C a , 7 a and if a only depend on the type of the incoming partons. 



The expression of Eq. ( 1791 coincides with the existing results for colourless final state NNLO 
cross sections, i.e., Drell-Yan process [H] and Higgs boson production via gluon fusion [191 EQ] 
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[j]. For the latter, when we consider the Higgs-gluon effective coupling in the large top mass (M t ) 
limit, the corresponding vertex factor C\ has the following perturbative expansion in powers of 
«s UlllII]: 

lias {®s\ 2 T2777 19, u 2 ( 67 1, a 2 R \] ml ^ . . 

Ca = 1 + + — -\ In + Nf( + - n -^-4 + Oial) . 81 

4 vr V Ti J I 288 16 M 2 t \ 96 3 Mf J J V s; V ; 

The "tree-level" cross section a is proportional to a\C 2 and then can be expanded in powers of 
as, being the lowest order term proportional to a|. The NLO and NNLO corrections are then 
of order a| and «g respectively. To obtain the correct result when using Eqs. (|3"7|) and (|79|) to 
compute Eq. (Ej), one has to include the terms in Eq. ljBTl) up to the desired precision in the strong 
coupling constant. 

The results we have presented can also be trivially expressed in a more general way, by undoing 
the formal (iPS^ integrals. Then one arrives at an identical expression for the differential cross 
section d 2 a/dQ 2 dPS^°\ where dPS^ depends on all the internal kinematical variables of the 
system F. In this way the result would be differential except on the hadronic activity which has 
been integrated out. 



5 Soft-virtual approximation in Mellin space 

The soft limit can be defined in a more natural way by working in Mellin (or iV-moment) space, 
where instead of distributions in z the dominant contributions are provided by continuous func- 
tions of the variable N. We consider the Mellin transform <Jn{Q 2 ) of the hadronic cross section 
ct(sh, Q 2 )- The iV-moments with respect to r = Q 2 /sh at fixed Q are thus defined as follows: 

a N (Q 2 ) = f drr N - 1 a(s H ,Q 2 ) . (82) 
Jo 

In iV-moment space, Eq. ([2]) takes a simple factorized form 

0"Af-l(Q 2 ) = O"0 ^2 fa,/hi,N{lA) fb/h 2 N(fJ 2 F) G a b, 7v( a S, Q 2 / Hr] Q 2 / lA) J (83) 



a,b 



where we have introduced the customary iV-moments of the parton distributions (f a /h,N) and of 
the hard coefficient function (G a &,/v) : 

fa/h,N^F) = dx l"" 1 fa/h(x,fip) , (84) 







1 

JV-1 



G abtN = / dz z™- 1 G ab (z) . (85) 





Once these iV-moments are known, the physical cross section in z-space can be obtained by Mellin 
inversion: 



/ ^ . V- [ Cmp+IQ ° dN fQ 2 y N+1 , 
ab Jc M p-ioo * m \s H ; 



x G o6iJV (a s ,Q7/4;<97/4) > (86) 



I Higher order soft (non- virtual) corrections to the T>i terms are obtained in P3J [351 [3j5] 
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where the constant Cmp that defines the integration contour in the iV-plane is on the right of all 
the possible singularities of the integrand. 



The evaluation of G a s, in the limit z — >■ 1 corresponds to the evaluation of the iV-moments 
G a a,N in the limit N — > oo. In the Mellin space soft- virtual approximation we drop all the terms 
that vanish when N — > oo and keep only constant and logarithmic (In N) contributions EE We 
introduce the notation SV-iV to indicate this approximation, while SV-z stands for the previous 
results obtained in the z space by keeping the most divergent terms when z — > 1. 

The (soft- virtual approximation to the) iV-moments G aa ,N can be obtained again as an expres- 
sion valid to all orders in e by using the following results: 



dz z 



N-l 



dz z 



N-l 



6(1 - z) 
1 



1 



T(N)Y{ 



-ae) 



T(N - ae) 



N° 



dz z N ' 1 V (z) 



(87) 
(88) 
(89) 



since all the contributions to the function G aS , have one of these three dependences on the variable 
z (the T>o(z) terms appear in the counterterms coming from mass factorization). Then, to obtain 
the SV-N approximation we just have to replace 5(1 — x) — > 1 and (1 — z)~ 1 ~ ae — > N ae T(—ae) in 
Eqs. ([2SD, (J3QD, (HO]), §HD, (USD and (J77D, and replace V (z) -> -ln(iV) -<y E in the Altarelli-Parisi 
splitting functions. 



We can also arrive at the SV-iV approximation from the expanded results of Eqs. (1371) and 
79\i . For this, we need the iV — > oo limit of the iV-moments of the distributions T>i(z), which 
can be obtained for instance from [I3l 42J. The results, keeping only the O(e ) terms, are the 
following: 



G 
G 



(i) 

aa,N 

(2) 
aa,N 



4 C a \n z (N) + 8 C alE ln(iV) + + O { - 



8C 2 Jn\N) + C a \n 3 (N) 



44 8 
32C alE + —C A --N f 



+ C a ln 2 (7V) 



- N f 



+ CJn(N) 



167T 2 



20 % 
9 3 



a 



C a [32 lE + 



+ 4- 
32 lE ir 2 



+ C A 

(!) n 
fin(O) 



134 2tt 2 



9 



+ 



44 7e \ 

3 ; 



+ ^(-14C 3 + ^ + 268 ^ 



56 , 40 7s , 8 7 | 



27 



9 



+ 



+ 8 7b 



(i) 

fin(O) 
O"0 



27 



9 



(90) 



44 7 | 



4 7E 7T 2 



(91) 



**Non-diagonal channels G a b,N with b a result in corrections which are at least O (-k) suppressed at the 
partonic level 
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where the coefficients and C^f are independent of A" and take the form 
4vr 2 . „ o o 



2 AW 

L-aa = -+4G q7e H , (92) 



00 



3 9 

, 55C 3 14 . 607 404 7s 134 7 2 44 7 | 67vr 2 2 7 ivr 2 37tt 
+ ^a^A 777 14 7sC3 + "777" H ^ 1 ^ 1 ^ h 



+ C a Nt I — — ^- — — — 



36 ' " 81 27 9 9 16 3 144 

5C 3 82 56 7s 20 7 2 8 7 i 5tt" 



18 81 27 9 9 
Hvr 3 „ o"£ } (o) /4vr 2 2 \ ^(o) 



6 cr V 3 / cr 



Again, the corresponding result for Higgs production using the ggH effective coupling has to 
be treated by incorporating the corrections described by Eq. (1ST]) . 



6 Phenomenological Results 

To evaluate the phenomenological accuracy of the SV-z and SV-A" approximations we compare 
them to the exact calculation for the processes of the kind of Eq. ([I]) that are available up to NNLO, 
i.e., Drell-Yan and Higgs boson production via gluon fusion. We compute the corresponding 
hadronic cross sections for proton-proton collisions at a center-of-mass energy A /s77 = 14 TeV. For 
simplicity, in the Drell-Yan process we only consider the photon channel, since the nature of the 
exchanged boson does not affect the impact of the QCD corrections we want to evaluate. 

To obtain the hadronic cross section we have to perform the convolution of the partonic result 
with the parton distribution functions. At each order, we use the corresponding MSTW2008 [13] 
parton distribution set and QCD coupling (one-loop «s at LO, two-loop as at NLO and three-loop 
a s at NNLO). 

We introduce now the notation used for the different contributions to the hadronic cross section. 
The up to NLO and NNLO calculations are denoted by <Jnlo(sh) and &nnlo{sh) respectively, 
while aLo( s H) is the LO cross section. The contribution cr^^(s^) is defined as the 0(a l s ) correction 
of the fixed order calculation ct wlo (sh), i.e., 

^lo(^) = ctloM + ■ ■ • + {^^M • (94) 

For the corresponding soft-virtual hadronic cross sections we use the same notation adding the 
index SV-z or SV-AT. The cross section ct^ lo (sh) contains the full NLO hard cross section. 

Since the soft-virtual cross section approximates only the dominant parton subprocess, to 
evaluate its accuracy we first have to compare it to the partial contributions for that partonic 
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Figure 1: The ratio between SV-N (solid lines) and SV-z (dashed lines) approximations and the 
exact result for Higgs boson production gg channel at the LHC (s/sh = 14 TeV) for the NLO 
(left) and NNLO (right) contributions. 



channel. In figured] we plot the quantities 

_ Vsv-zM K {i) ^sVjv(fg) rQ KN 



for i — 1, 2. In the left-hand side of figured] we show K^j_ z and K^J_ N for Higgs boson production 

(2) 

as a function of the Higgs mass Mr- = Q- The corresponding comparison at the next order K^_ z 

(2) 

and K§y_ N is shown in the right-hand side of the same figure. We present the equivalent plots 
for the Drell-Yan process in figure [21 as a function of the lepton pair mass Q. In all figures the 
central curves are obtained by fixing hr = = Q, and the bands by varying simultaneously the 
renormalization and factorization scales to [ar = [af = Q/2 and lir = lif = 2Q- 

We can see from the plots that in all the cases the iV-space approximation is in very good 
agreement with the exact NLO and NNLO calculations, being much more accurate than the z- 
space one, as expected since the soft- virtual approximation is formally defined is a more systematic 
way in Mellin space. For Higgs boson production the SV-N approximation reproduces about 90% 
of the exact result at each order, while for the Drell-Yan process the agreement reaches the level of 
80%. This difference is expected because the gluon distribution function grows faster than quark 
distributions for small fractions of x, enhancing the threshold contribution to the cross section for 
gluon fusion processes. 

Our next step is to evaluate if the inclusion of the N*LO soft-virtual corrections results in an 
improvement over the N i_1 LO calculation. We plot the ratio between 0"^Y£o( s #) an< ^ °"n i lo( s #)' 
and compare it with the ratio between cr N i-i LO (sH) and 0"n^lo( s #)> f° r £ = 1> 2. This is shown for 
Higgs boson production in figure [3] and for the Drell-Yan process in figure H] 
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Figure 2: The ratio between SV-N (solid lines) and SV-z (dashed lines) approximations and the 
exact result for the Drell-Yan qq channel at the LHC (y/sH = 14 TeV) for the NLO (left) and 
NNLO (right) contributions. 
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Figure 3: The ratio between <7^Y L o (solid lines) and cr N i-i LO (dashed lines) and cr N i LO for Higgs 
boson production at the LHC (\/sh = 14 TeV), for i — 1 (left) and i = 2 (right). 
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Figure 4: The ratio between o"^Y L q (solid lines) and cr N i-i LO (dashed lines) and cr N ; LO for the 
Drell-Yan process at the LHC (s/sh = 14 TeV), for % — 1 (left) and i — 2 (right). 

For Higgs boson production we see that the ratio between the SV-iV approximation and the 
full calculation is very close to 1 for both NLO and NNLO [13]. We can also see that, due to the 
large perturbative corrections, the N i_1 LO calculation is far from the N*LO one (for i = 1, 2). 
If the full NNLO result were not available, it would be clearly convenient to include the SV-iV 
approximation contribution as an attempt to improve the accuracy of the calculation. Given 
the fact that one can reach this conclusion a priori, simply based on the relevance on soft-gluon 
emission in gluon initiated processes, it seems very reasonable to conjecture that the soft-virtual 
approximation might be equally accurate for any gluon-gluon fusion dominated reaction. 

The situation is different for the Drell-Yan process, mainly for two reasons. In the first place, 
the total cross section has more sizeable contributions from other partonic subprocesses besides 
the quark- antiquark one, mainly from the quark-gluon channel which is not included in the SV 
approximation. In the second place, the perturbative corrections are small compared to those for 
Higgs boson production, and then the N l_1 LO calculation is closer to the N*LO one. However, we 
can see that the NLO SV-iV approximation is far more accurate than the LO result, and the NNLO 
SV-iV one is as accurate as the NLO result, improving with the invariant mass. Therefore, one 
might expect even more accurate results for quark initiated process where heavier states (closer 
to the production threshold for the same collider energy) are produced in the final state, as in 
double gauge boson production. 



7 Threshold Resummation 

Given that soft-virtual terms provide the bulk of the corrections for the processes under study, it 
is possible to improve over the state of the art fixed-order predictions by performing soft-gluon 
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resummation. In this section we consider the all-order perturbative summation of enhanced thresh- 
old (soft and virtual) contributions to the partonic cross section. We refrain from analyzing the 
phenomenological impact of soft-gluon resummation and, instead, concentrate on the extraction 
of a universal expression for the coefficients needed to achieve next-to-next-to-leading logarithmic 
(NNLL) accuracy. 

The formalism to systematically perform soft-gluon resummation for hadronic processes, in 
which a colourless massive system F is produced by qq annihilation or gg fusion, was set up in 
Refs. [II |45l|46]. 

The resummation of soft-gluon effects is achieved by organizing the partonic coefficient function 
in Mellin space as 

gS(«s(^), MZ/vk K/ti) = 3*(«s(a&), <//4) 

• A N (a s (v 2 R ), M 2 H /n 2 R ; M 2 H /&) + 0(1/N) , (96) 

The large logarithmic corrections (that appear as \n 2n ~ k iV in Mellin space) are exponentiated 
in the Sudakov radiative factor A^v, which depends only on the dynamics of soft gluon emission 
from the initial state partons. It can be expanded as 

A N as (/£), In N; -f, -f )= In N g^(Qa s (^ R ) In N) + #(Qa s (/4) In N, M 2 H /fi 2 R ; M 2 H /fi 2 F ) 

+ <*s(A) 9?HQ<xs(jJ?r) ln ^ M h/^ M 2 H /fi 2 F ) 
+ E [«s(^)] n_2 ^ n) (Q«s(/i|) IniV, M 2 H /fi 2 R ; M 2 H /^ 2 F ) . (97) 

n=4 



The function ln N ga resums all the leading logarithmic (LL) contributions «g ln n+1 N, g^ 
contains the next-to-leading logarithmic (NLL) terms ag ln n N, a^ga collects the next-to-next- 
to-leading logarithmic (NNLL) terms ag +1 ln n N, and so forth. All the perturbative coefficients 
required to construct the gcP,g^ and g^ functions are known and only depend on the nature of 
the initiating partons. Their explicit expression can be found, for instance, in Ref . [T3l I4"7] . 

On the other hand, the function C aa {cts) contains all the contributions that are constant in 
the large- iV limit. They are produced by the hard virtual contributions and non- logarithmic soft 
corrections, and can be computed as a power series expansion in as'- 

CU*s(»l),M 2 H /^M 2 H /» 2 F ) = l + J2 CS{M 2 H /^M 2 H /^) . (98) 

n=l ^ 71 ' 



The Cai coefficient required to perform the resummation up to N l LL can be obtained from the 
corresponding fixed order computation to N l LO accuracy. No general expression was known up 
to now for the hard coefficient, and its obtention demanded an individual computation for each 
process. We show here that, given the results for the soft- virtual approximation in Mellin space 
presented in the previous section, it is possible to obtain the universal expressions for the C^a 
coefficients up to NNLL accuracy. 
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A direct comparison of the expansion of Eq. (l9"oT) up to order «g to the expression presented 
in Eqs. (1901) and fl9T|) shows a complete agreement for the logarithmically enhanced terms (which 
are fully predicted by the resummed expression), and allows the extraction of the C^l coefficients 
from the N— independent terms. As a matter of fact, since it is customary to define the functions 
ga such that (Octs In N) = when as = 0, the Sudakov radiative factor can only produce 
logarithmically enhanced terms (i.e., no constant terms) when expanded to any fixed order in the 
strong coupling constant and, therefore, the C^a are exactly given by the N— independent terms 
of the soft-virtual coefficients G^l N at each order in perturbation theory. Up to NNLL accuracy 
they are expressed by Eqs. (l92l and ([93]) as 

Mi) _ r (i) 

^ aa ^ aa 5 

cS = cS ■ (99) 

The coefficients in Eq. (1991) depend, as expected, on the virtual corrections to the corresponding 
scattering amplitudes and on the colour factors of the initial state partons. We have explicitly 
checked that the known coefficient for Higgs production [13] agrees with the result from Eq. (l99|) . 

The new universal expression in Eq.( l99l) plus the present knowledge of the Sudakov radia- 
tive factor allows us to perform the resummation of soft-gluon emission relevant at the partonic 
threshold up to NNLL accuracy for any process of the kind h\ + h 2 — > F + X. 



8 Conclusions 



In this paper we have computed the NNLO soft and virtual QCD corrections for the partonic cross 
section of colourless-final state processes in hadronic collisions. We presented a universal expres- 
sion for the corresponding cross section, whose only dependence on the process enters through the 
(finite part of) one- and two-loop amplitudes. 

With the recent calculation of the three- loop quark and gluon form-factors [HI 09[ [50], it 
should be possible to attempt for an evaluation of the dominant soft-virtual corrections at N 3 LO 
for a number of interesting processes (as Drell-Yan and Higgs boson production) following the 
approach described in this paper. As a first step in this direction, we have provided the necessary 
ingredients for ultraviolet and infrared factorization by presenting explicit results for the cross 
section valid to all orders in e. 

We evaluated the accuracy of the soft-virtual approximation for known processes as Drell-Yan 
and Higgs boson production in hadronic colliders. We conclude that the approximation is excellent 
for gluon-gluon fusion initiated processes and quite accurate for quark-antiquark initiated ones 
for high invariant masses. For processes for which it is still not possible to compute the full 
NNLO corrections, counting with the corresponding soft-virtual approximation results in a clear 
improvement over the accuracy of the available calculation. 

Finally, profiting from the soft-virtual calculation, we provided a universal expression for the 
coefficient needed to perform threshold resummation up to NNLL accuracy. 

tt After properly taking into account the corrections to the effective ggH vertex 



22 



Acknowledgements 

This work was supported in part by UBACYT, CONICET, ANPCyT and the Research Executive 
Agency (REA) of the European Union under the Grant Agreement number PITN-GA-2010-264564 
(LHCPhenoNet). 

References 

[1] R. Hamberg, W. L. van Neerven and T. Matsuura, Nucl. Phys. B 359, 343 (1991) [Erratum- 
ibid. B 644, 403 (2002)]. 

[2] R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett. 88, 201801 (2002) |hep-ph/020i206] . 

[3] C. Anastasiou and K. Melnikov, Nucl. Phys. B 646, 220 (2002) |hep-ph/0207004] . 

[4] V . Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B 665 (2003) 325 
[hep-ph/0302135] . 

[5] S. Dawson, Nucl. Phys. B 359 (1991) 283. 

[6] A. Djouadi, M. Spira and P. M. Zerwas, Phys. Lett. B 264 (1991) 440; M. Spira, A. Djouadi, 
D. Graudenz and P. M. Zerwas, Nucl. Phys. B 453 (1995) 17. 

[7] M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, Nucl. Phys. B 453 (1995) 17. 

[8] D. de Florian and M. Grazzini. larXiv:1206.413"3l [hep-ph], 

[9] J. Baglio and A. Djouadi, JHEP 1103 (2011) 055. 

[10] C. Anastasiou, S. Buehler, F. Herzog and A. Lazopoulos, JHEP 1204, 004 (2012). 

[11] D. de Florian and M. Grazzini, Phys. Lett. B 674, 291 (2009). 

[12] C. Anastasiou, R. Boughezal and F. Petriello, JHEP 0904 (2009) 003. 

[13] S. Catani, D. de Florian, M. Grazzini and P. Nason, JHEP 0307 (2003) 028 [hep-ph/03062lT] . 

[14] S. Moch and A. Vogt, Phys. Lett. B 631 (2005) 48 |hep-ph/0508265] . 

[15] E. Laenen and L. Magnea, Phys. Lett. B 632 (2006) 270 |hep-ph/0508~284] . 

[16] V. Ravindran, Nucl. Phys. B 752 (2006) 173. 

[17] M. Kramer, E. Laenen and M. Spira, Nucl. Phys. B 511 (1998) 523. 

[18] T. Matsuura, S. C. van der Marck and W. L. van Neerven, Nucl. Phys. B 319, 570 (1989). 

[19] S. Catani, D. de Florian and M. Grazzini, JHEP 0105, 025 (2001) |hep-ph/0102"227] . 

[20] R. V. Harlander and W. B. Kilgore, Phys. Rev. D 64 (2001) 013015 |hep-ph/0102'24l] . 



23 



[21] R. K. Ellis, M. A. Furman, H. E. Haber and I. Hinchliffe, Nucl. Phys. B 173 (1980) 397. 
[22] A. Bassetto, M. Ciafaloni and G. Marchesini, Phys. Rept. 100 (1983) 201. 
[23] W. T. Giele and E. W. N. Glover, Phys. Rev. D 46, 1980 (1992). 

[24] S. Catani and M.H. Seymour, Phys. Lett. B 378, 287 (1996), Nucl. Phys. B 485, 291 (1997) 
(E ibid. B 510, 503 (1998)). 

[25] Z. Kunszt, A. Signer and Z. Trocsanyi, Nucl. Phys. B 420, 550 (1994) |hep-ph/940i~2~9~4] . 

[26] S. Catani, Phys. Lett. B 427, 161 (1998) [hep-ph/9802439] . 

[27] R. K. Ellis, W. J. Stirling and B. R. Webber, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 
8 (1996) 1. 

[28] S. Catani and M. Grazzini, Nucl. Phys. B 591, 435 (2000) [hep-ph/0007142] . 

[29] Z. Bern and G. Chalmers, Nucl. Phys. B 447 (1995) 465 |hep-ph/9503'236] . 

[30] Z. Bern, V. Del Duca and C. R. Schmidt, Phys. Lett. B 445 (1998) 168 lhep-ph/9810~40~9~] . 

[31] Z. Bern, V. De l Duca, W. B. Kilgore and C. R. Schmidt, Phys. Rev. D 60 (1999) 116001 
[hep-ph/9903516] . 

[32] S . Catani, D. de Florian, G. Rodrigo and W. Vogelsang, Phys. Rev. Lett. 93 (2004) 152003 
[hep-ph/0404240] . 

[33] S. Catani and M. Grazzini, Nucl. Phys. B 570, 287 (2000) |hep-ph/9908523] . 

[34] F. A. Berends and W. T. Giele, Nucl. Phys. B 313 (1989) 595. 

[35] C. Anastasiou, S. Buehler, C. Duhr and F. Herzog. TarXiv: 1208.3 130 [hep-ph]. 

[36] S . M. Aybat, L. J. Dixon and G. F. Sterman Phys. Rev. D 74 (2006) 074004 
[hep-ph/06073 09|; Phys. Rev. Lett. 97 (2006) 072001 |hep-ph/ 06 06254] . 

[37] Z. Bern, L. J. Dixon and D. A. Kosower, JHEP 0408, 012 (2004) |hep-ph/0404"293| . 

[38] V. Ravindran, Nucl. Phys. B 746 (2006) 58 |hep-ph/05 12249] . 

[39] V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B 767 (2007) 100 
[hep-ph/0608308] . 

[40] K. G. Chetyrk in, B. A. Kniehl and M. Steinhauser, Phys. Rev. Lett. 79, 353 (1997) 
[hep-ph/9705240] . 

[41] M. Kramer, E. Laenen and M. Spira, Nucl. Phys. B 511, 523 (1998) |hep-ph/ 96 11272] . 
[42] S. Catani and L. Trentadue, Nucl. Phys. B 327, 323 (1989). 

[43] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 
[arXiv:0901.0002l [hep-ph] ] . 



24 



[44] G. Sterman, Nucl. Phys. B 281 (1987) 310. 
[45] S. Catani and L. Trentadue, Nucl. Phys. B 327 (1989) 323. 
[46] S. Catani and L. Trentadue, Nucl. Phys. B 353 (1991) 183. 
[47] A. Vogt, Phys. Lett. B 497 (2001) 228 [hep-ph/0010146) . 

[48] T. Gehrmann, E. W. N. Glover, T. Huber, N. Ikizlerli and C. Studerus, JHEP 1006 (2010) 
094 jarXiv: 1004.36531 [hep-ph]]. 

[49] R. N. Lee, A. V. Smirnov and V. A. Smirnov, JHEP 1004 (2010) 020 [ arXiv:1001.2887l 
[hep-ph]]. 

[50] P. A. Baikov, K. G. Chetyrkin, A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. 
Rev. Lett. 102 (2009) 212002 jarXiv:0902.35l9l [hep-ph]]. 



25 



